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Abstract 

I construct a two-dimensional lattice on which the inhomogeneous site percolation threshold is 
exactly calculable and use this result to find two more lattices on which the site thresholds can 
be determined. The primary lattice studied here, the "martini lattice", is a hexagonal lattice 
with every second site transformed into a triangle. The site threshold of this lattice is found to 
be 0.764826..., while the others have 0.618034... and l/\/2. This last solution suggests a possible 
approach to establishing the bound for the hexagonal site threshold, pc < l/\/2- To derive these 
results, I solve a correlated bond problem on the hexagonal lattice by use of the star-triangle 
transformation and then, by a particular choice of correlations, solve the site problem on the 
martini lattice. 
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FIG. 1: The martini lattice, drawn to emphasize its origin as a hexagonal lattice with every second 
site transformed into a triangle. 



I. INTRODUCTION 

The star-triangle transformation has been a useful tool in 2-d percolation theory 
was originally used by Sykes and Essam ^ , along with planar duality, to find the thresholds 
for inhomogeneous bond percolation on the triangular and hexagonal lattices, and was later 
adapted by Wierman ^ to solve the bond problem on the bowtie lattice and its dual. In 
fact, all exact results in 2-d are either derived directly from duality or matching properties, 
or rely on the star-triangle transformation in some way. In the present work, this method 
will be extended to problems with some limited correlation structure. 

The primary lattice studied here, which I call the "martini lattice" due to the shape of 
the basic cell, is the one shown in figures Q and El Each site has 3 nearest neighbors, but 
the lattice is non-uniform because some sites are (3, 9^) while others are (9^) in the notation 

n 

of Griinbaum and Shephard p]. This lattice is mentioned on page 186 of that book as an 
example of what they call a 2-homeohedral tiling of valence 3. 



II. STAR-TRIANGLE TRANSFORMATION 



The star-triangle transformation exploits the fact that if the bonds of a unit cell of the 
triangular lattice (T) are replaced by a corresponding star, as illustrated in figure El the 
result is the hexagonal lattice (H), which is the triangular lattice's dual (figure |31). The bond 
percolation threshold of a lattice, L, and its dual, L'^, are related, in 2-d, by the well-known 
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FIG. 2: The martini lattice, drawn to emphasize the martinis. 




A ^ B 

FIG. 3: The star-triangle transformation, p, r, s denote the probabilities of their corresponding 
bonds on the triangle and star, and A, B, C label the sites. 

formula 

This means that the appearance of the infinite open cluster on L coincides with the dis- 
appearance of the infinite closed cluster on L'^. The star-triangle transformation leads to 
another relationship between the critical probabilities of H and T besides Pc(T) = 1— Pc(H), 
which allows both to be determined exactly. The method even works for the inhomogeneous 
case, where the probabilities of each bond being open on the base triangle are different, 
resulting in a critical surface rather than a critical point. 

The argument proceeds as follows. Consider bond percolation on the triangle and super- 
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imposed star shown in figure El The probabihties p,s,r refer to the probabilities that their 
corresponding bonds are open on either the star or triangle. We can ask several questions 
about the connectedness of the sites A, B, C . For example, what is the probability that A 
is connected to both B and C, an event we will denote P{A ^ B, A ^ C), through open 
bonds on the triangle? This is easily found to be 

P{A ^ B,A^ C) =ps+ pr{l -s) + sr(l - p) (2) 

Next, we want the probability that A,B,C are connected through closed bonds on the star. 
We denote this event Q*{A B, A ^ C). Q will hereafter denote the probability of events 
that happen in closed bonds and * will indicate that the event happens by traversing the 
star rather than the triangle. Since this event only happens if all three bonds are closed, we 
have 

Q*^A^B,A^C) = il- p)il - r)(l - s) (3) 

If we now consider the entire lattice, we can see that the replacement of triangles by stars 
turns the triangular lattice into the hexagonal lattice, i.e. the dual (figure |3)). This implies 
that the condition 

P{A^ B,A^C) =Q*{A^ B,A^C) (4) 

defines our critical surface. The connectivity of open bonds on T is exactly the same as that 
of the closed bonds on H when this condition is satisfied. This means that if there is an 
infinite open cluster on T, then there is an infinite closed cluster on H. This leads, by our 
previous discussion of duality, to the conclusion that there is neither an infinite open cluster 
on T nor an infinite closed cluster on H, i.e. we are on the critical surface. Substituting our 
results into (jlj and simplifying, we are led to the final result: 

psr — p — s — r + l = (5) 

Of course, setting s = r = p leads to the critical point for bond percolation on the triangular 
lattice, pc = 2sin7r/18. We can also check that the same result is obtained by studying 
different connectivities of the sites. For example, the probability that A connects B, but 
not C, denoted P{A ^ B,A ^ C), is given by: 

P{A^ B,A^C) =p{l-p)'^ (6) 
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Also, 

Q*{A^B,A^C)=pil-py (7) 

so that P{A —>■ B, A ^ C) = Q*{A ^ B, A ^ C) for all p, so this does not provide any 
constraint. If we calculate P{A ^ B, A ^ C) = Q* {A ^ B , A ^ C) , we recover (0) as we 
should. 

III. CORRELATED BOND PERCOLATION ON THE TRIANGULAR LATTICE 

The method of Sykes and Essam can be extended to the case where the bonds of the 
triangle are not independent. Our critical surface will now appear as a constraint between 
the 1, 2, and 3-point joint probabilities. It is important to note that although all bonds in 
a triangle are correlated, there are no correlations between neighboring triangles, so a given 
bond is only correlated to two of its neighbors. The dual lattice is constructed in the same 
way as in the uncorrelated case, with each bond in the dual inheriting the probabilities and 
correlations of the original lattice. Labelling the bonds v, h, I as shown in figure we will 
deal with the quantities P{h,v,l), P{h,l), P{v,h), P{v,l), P{v), P{h), P{1), which are the 
set of 1, 2, and 3-point joint probabilities of the indicated bonds being open. Probabilities 
of bonds being closed will be denoted with a bar over the bond name, e.g. P{v). We can 
now repeat the procedure outlined above but with our joint probabilities: 

P{A -^B,A^C) = P{v, h) + P{v, I, h) + P(/i, /, v) (8) 

and 

Q*{A^ B,A^ C) = P{v, I, h) (9) 
Equating these gives our critical surface: 

P{v, h) + P(v, /, h) + P{h, /, v) - P{v, lh)=0 (10) 

There are many equivalent ways this can be expressed. For example, if we use the condition 

P{A ^ B,A^C) =Q*{A^ B,A^C) 
we obtain the more compact 

P{v) + P{v,h,l) - P{hJ) = (11) 
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FIG. 4: The star-triangle transformation. Replacing each triangle with a dashed star transforms 
the triangular lattice into the hexagonal lattice. 




FIG. 5: Labels used in treating the correlated triangle. A, B, C label sites and w, /i, I label bonds 
on the triangle and the corresponding bonds on the star - note that these are now names not 
probabilities. 

Although they look dissimilar, equations (jlU|) and (jllj) are in fact the same constraint and 
it is a simple matter to relate them using identities. 

To compare with our earlier results, it is easy to see that setting P(/i, /) = (! — r)(l — s), 
P{v, h, Z) = (1 — p)rs, P{v) = p leads to 
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FIG. 6: The assignment of probabilities to the sites on the martini lattice. 

IV. SITE-TO-BOND TRANSFORMATION 

If we consider a realization of site percolation on a given lattice, we can transform it into 
a bond process by declaring a bond to be open if both its bounding sites are occupied, and 
closed otherwise. By doing this we introduce correlations between neighboring bonds; the 
probability that a given bond is open is p^, but the probability that a bond is open given 
that one of its neighbors is open is p. Furthermore, there are also three-point correlations; 
the probability that a bond is open is 1 if two of its neighbors on opposite ends of the bond 
are open. It is clear that the existence or lack of an infinite open cluster are properties that 
will be shared by both the site and transformed bond problems. If we now consider sites on 
a triangle, we can use these rules to derive joint probabilities for the bonds, which we can 
then use in the criticality condition (fTTj). It is easy to see that if the sites are occupied with 
probability p: 

P{v) = p^ (12) 
P{hJ) = {I - pf + "ipil - pf + p\l - p) (13) 
P{h,l,v) = (14) 

If we use these in we do not expect to have solved the site problem on the triangular 
lattice. The critical surface is not appropriate to that problem because we have not included 
correlations between triangles. I suggest that the threshold we will discover is that of the site 
problem on the Kagome lattice (figureHj). To see this, consider the triangles outlined in figure 
|7| Clearly they are not correlated to each other if we use the site-to-bond transformation. 
But percolation of bonds on these triangles implies percolation of the lattice since, due to 
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3-point correlations, the bonds on the separating triangles will be open with probability 1 if 
two bonds on either side of it are open. Plugging expressions ()12|) - |T^ into (|TT|) we obtain 
the polynomial 

1~3pI+pI = (15) 

with solution pc = 1 — 2sin7r/18 which is indeed the critical threshold of the Kagome 
lattice. However, since this is just the covering lattice of H, its threshold has long been 
known through more elementary means. 

We can obtain our new results by considering percolation on the star. The critical surface 
in this case will be given by the complement of (fTTj) : 

P{v) + P{v,hJ) - P{hJ) = (16) 

In fact, we will consider the inhomogeneous site problem, and assign probabilities p,r,s,t, to 
the sites as shown in figure El 

P[v) = 1-st (17) 
P{h,l) = prs (18) 
P{v,h,l) = st{l-p){l-r) (19) 

This leads to the critical surface 

1 — rst — prs — pst + stpr = (20) 

which is the central result of this work. But to what lattice does it correspond? In the 
previous example, where we obtained the Kagome lattice, we needed to insert extra triangles 
to separate our correlated triangles. Inserting these separating triangles in between the stars, 
we obtain the martini lattice shown in figure |21 The critical threshold for site percolation is 
obtained by setting r = s = t = p: 

l-3pl+pt = (21) 

which has solution on [0, 1]: 

Pc = 0.764826... (22) 

We can obtain further results by making different choices for the probabilities. For example, 
by setting s = l,r = t=p, we expect the Kagome lattice to reappear, since s = 1 turns 
the star back into a triangle. Plugging these into (p?)) we indeed get (fT3)). Other choices are 
possible that lead to a variety of results. 
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FIG. 7: The Kagome lattice. The circled triangles are the ones on which we apply the site-to-bond 
transformation. 




FIG. 8: The lattice obtained from the martini lattice by setting t = r = 1. On this lattice, 
pf^ = {^/5-l)/2 = 0.618034... and = 1/2. 

A. t = r = 1 

The corresponding lattice is the one shown in figure |H1 and it resembles a stack of houses 
or a neighborhood. It is non-uniform, with some sites (3^, 5^) with 5 nearest neighbors and 
others (3, 5^) with 3 nearest neighbors, and falls somewhere between the hexagonal lattice 
and the (3^,4^) lattice - almost exactly between as we will see. The (3^,5^) sites have 
probability s and the (3, 5^) sites p. The critical locus is 

l-s-ps = (23) 

Setting s = p leads to the critical threshold 

1-Pc-pI = (24) 
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or = (v^ - l)/2 = 0.618034... 

This lattice has some interesting properties that are worth mentioning. For one, it is 
self-dual, meaning we can immediately locate its bond threshold at 1/2, as well as the site 
threshold of its covering lattice. Also, it seems to be an intermediate step in transforming the 
hexagonal lattice into the (3^,4^) lattice. If we remove the horizontal bonds, the "floors" of 
the houses, we get the hexagonal lattice. If instead we symmetrically add another horizontal 
)ond, a "crossbeam" under the roof, we get the (3^, 4^) lattice. It has been found numerically 
^ that Pc(33,42) ^ 0.550 and pdU) ^ .697 and their average is around 0.624, which is 
similar to the result for the present lattice. Apparently, adding the floor has roughly the 
same effect on the threshold as adding the crossbeam. 

B. s = t = l 

Setting s = t = 1 and r = p, we get the covering lattice of the square bond problem, 
leading to Pc = 1/2, which is a very roundabout way of solving that problem. 

C. r = 1 

Again, a star is turned into a triangle, but a different one from that which produced the 
Kagome lattice earlier. The lattice that results here is shown in figure IHl The critical surface 
is 

l-st-ps = (25) 

Setting s = t = p yields 

l-2pl = (26) 

which means pc = I/V2 = 0.707.... This is an interesting result for several reasons. For 
one, some sites have 3 nearest neighbors while others have 4. Thus, if we were to make a 
guess strictly on the basis of nearest neighbors, we might be led to believe that this lattice's 
site threshold was smaller than that of the hexagonal lattice, where every site has 3 nearest 
neighbors. The contrary is true however, as we know from numerical results j^. Interesting 
Pc = was once conjectured to be the exact site threshold for the hexagonal lattice 

but was shortly thereafter judged unlikely from numerical considerations 
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FIG. 9: The lattice obtained by setting r = 1. The site threshold is pc = \j\f2. 




FIG. 10: The transformation that takes the hexagonal lattice to the one shown in figure |H1 The 
site at the bottom of every hexagon is divided into two sites, and a bond is inserted between them. 

This result suggests an approach to finding an upper bound for p^'**'(H). Consider the 
following procedure for producing this lattice. Starting with the usual hexagonal lattice, take 
the site at the bottom of each hexagon, split it into two sites and connect them together 
with a bond, as shown in figure^) All that is required is to show that this procedure always 
increases the critical probability and we will have shown that 

pf ^(H) < i= (27) 

Of course, even better would be if the shift in the critical point could be exactly determined. 
It remains to be seen whether either of these suggestions is actually feasible. 

V. CONCLUDING REMARKS 

We have found three lattices whose site percolation thresholds can be calculated exactly, 
and shown that one of these solutions might lead to an improved upper bound for the 
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hexagonal lattice. 

It should be noted that a special case of has previously appeared in the literature, 
ihough in a slightly different context. In considering a mixed site/bond problem, Kondor 
3] derived an expression that matches (j2(jp for r = t = p. The situation he considered was 
bond percolation on the triangular lattice in which there is an "activating" site in the middle 
of every second triangle that turns all its surrounding bonds off if it is unoccupied. This 
problem is evidently isomorphic to site percolation on the martini lattice. 

Also, the value 1/V2 for the lattice in figure IHl is actually somewhat puzzling. Usually, 
simple values like this indicate an elementary tranformation from a known lattice with 
Pc = 1/2. If such a transformation exists in this case, it does not appear to be obvious. 
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